The 1-d stochastic wave equation 
driven by a fractional Brownian motion 

'^ ; Lluis Quer-Sardanyons * 



O 
O 



O 
-(— > 

s 

X 



Project OMEGA, INRIA Lorraine 



pq ' Institut Elie Cartan, Campus Scientifique 

5_^ ! BP 239 - 545O6 Vandoeuvre-les- Nancy, France 

Qh' quer@iecn.u-nancy.fr 

'• Samy Tindel 

Institut Elie Cartan, Universite de Nancy 1 
BP 239 ~ 545O6 Vandceuvre-les- Nancy, France 
^ ' tindel@iecn.u-nancy.fr 

t^ : 8th February 2008 



^ , Abstract 

'NT 

f~^ I In this paper, we develop a Young integration theory in dimension 2 which will 

.^ ■ allow us to solve a non-linear one dimensional wave equation driven by an arbitrary 

^^ , signal whose rectangular increments satisfy some Holder regularity conditions, for 

some Holder exponent greater than 1/2. This result will be applied to the infinite 

dimensional fractional Brownian motion. 
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1 Introduction 

During the last past years, some spectacular advances have been made in order to define 
and solve some stochastic PDEs driven by a Brownian noise with a very general spatial 
covariance structure. This has been achieved for instance in the infinite dimensional 
setting in j^H 1201 ^^^ both heat and wave equations, while the multiparametric setting 
has been treated first mainly for the 1-d wave equation in jSHl and [2], and for the 2- 
dimensional case in jH] and ^\. The Walsh setting for SPDEs (J2S]) has been generalized 
then in |1], leading to a quite complete picture for the stochastic heat equation, and 
also for the wave equation up to dimension 3, allowing some deeper study in J22ll23ini- 
Notice also that an effort has been made in order to cover the case of Levy noises e.g. in 

1211 mini IE!- 



On the other hand, since the recent introduction of fractional integrals [29] and rough 
paths (see [IHIIISI) techniques in probability theory, it has been clear that the pathwise 
method could be a good way to extend the notion of stochastic differential equations to 
a wide variety of noises. However, while this strategy has been explored thoroughly in 
the case of ordinary differential equations (see e.g. [21101), the case of stochastic PDEs 
is still widely open. Indeed, if the case of linear heat |2H] or wave PP equations driven 
by fractional noises has been considered, sometimes leading to optimal results, only some 
very partial results are available in the case of non-linear equations: let us mention for 
instance ^Hl for stochastic evolution equations driven by a fractional noise, with a very 
regular space covariance. In this context, we have started to build in pQ3j an extension 
of Young integrals to the evolution setting, which allowed us to solve the stochastic heat 
equation driven by an infinite-dimensional fractional Brownian noise with a non-trivial 
spatial covariance structure. Observe that, while limited to the Young case, it seems that 
the methods introduced in the latter reference can be extended to the rough path case 
(see m)- 



The aim of the current paper is to show that the approach initiated in ^\ can be 
extended to hyperbolic PDEs, and we will illustrate this fact by considering the stochastic 
wave equation in M, of the form 



^'^) - 1^(^'^) = ^{y{s,t))X{s,t), for (s,t) e [0,T] X 



with initial conditions given by 

dY 

r (0, t) = -—(0, t) = 0, for t e M. 

OS 

In equation (^, a stands for a smooth function from M to M, and X represents the noise 
which drives the equation. As usual in the SPDE theory, (0) is understood in the mild 
sense, which can be specified as follows: we will say that y is a solution to (Q) if, for any 

(s, t) e [0, T] X M, we have 

Y{s,t)= f f a{Y{u,v))X{du,dv), (2) 

J JC{s,t) 



where C{s,t) denotes the hght cone with vertex (s,t), that is the triangle dehmited by 
the points (s,t), (0, t + s) and (0,t — s), and where the integral defining equation (J2)) is 
understood in the Young sense. Then, for this latter equation, we will give some existence 
and uniqueness results for a general class of noises whose rectangular increments are 
Holder continuous with Holder exponent > ^ (see Theorem 13 .41 for a precise statement). 

One particular case of interest for us will be the infinite-dimensional fractional Brow- 
nian motion, which can be defined in the following way: on a given complete probability 
space (r2,jF, P), let X be a centered Gaussian family {X(0);0 G ViM."^)} indexed by the 
set of test functions I?(]R^), with a covariance function given by 

E iXicfyXii))) = ch ! dudv\u~v\^^-^ f dxdy\x-y\-''(p(u,xU(v,y), (3) 

J[0,T]2 JR2 

where H > 1/2, ch = H{2H—1), and u G (0, 1). With a slight abuse of notation, set then 
X{s,t) = X(l[o,t]x[o,s])- Let also TP'^^D) be the space of functions defined on D d M? 
having a Holder regularity of order 7 in time and 7 in space (see Definition 12.11 for the 
precise requirements). Eventually, let TZ^il be the -45° degree rotation in the plane. In 
this context, our existence and uniqueness result will be the following: 

Theorem 1.1. Let X he a fractional Brownian noise defined by (0). Suppose that a is a 
bounded function in C^{M.) with bounded derivatives of any order and let f],fi & (0, 1) be 
such that ri + fj<l + H— ^. Then, Equation (OP admits a unique solution Y such that 
Y o TZ-iL G TiP'^{R), for any rectangle R around the origin. 

As mentioned above, this paper can be seen as an extension of ^0], but the method- 
ology used here is quite different from the evolution type considerations contained in this 
latter reference. Indeed, it seemed easier, in the case of the wave equation, to consider 
the problem at hand in the multiparametric setting. This has lead us to the following 
global strategy: 

1. Construct first a general Young integral on rectangles i? C M^ whose sides are par- 
allel to the axes. This integral will have the form /^ fdg., for two Holder continuous 
functions f , g : R —>■ M, with large enough Holder indexes. Notice that our con- 
struction is inspired by [Oj, but it is expressed here directly in terms of convergence 
of Riemann sums, while [5] uses a 2-d analog of the A- map defined in [S] . 

2. Extend this Young integral in order to cover the case of a domain R which is 
a triangle with two sides parallel to the axes. This is done in a straightforward 
manner, by writing the triangle as a countable union of rectangles. 

3. Rotate the wave equation in order to deal with an ordinary differential equation in M^ 
involving the previous triangular domains. Once our Young integral is constructed, 
the existence and uniqueness result will be obtained by an extension of the usual 
fixed point argument for differential equations. It is worth noticing here that our 
computations for this step will be quite delicate, in spite of having chosen a very 
regular coefficient a. Indeed, though a G C^(M), we will see that its interpretation 



as a map from TP'^ into itself does not enjoy the properties one usually assumes for 
the resolution of Young equations: in fact, it is only locally Lipschitz with quadratic 
growth, a fact which will add some technical difficulties to our analysis. 

4. In order to handle the case of the fractional Brownian noise, one has to show that 
the rotation of this noise still satisfies the Holder regularity conditions allowing the 
definition of a Young integral. This can be done in our case, thanks to some almost 
explicit and cumbersome calculations. 

This strategy will be made more explicit in the remainder of the paper, but let us mention 
at this point that, to our knowledge. Theorem 11.11 is the first existence and uniqueness 
result for a non-linear wave equation driven by a general kind of noise, and in particular 
by an infinite dimensional fractional Brownian motion. We hope to extend this approach 
to a more irregular noise in a subsequent publication. Let us also mention that some of 
our techniques can be related to those developped in J2ZI for numerical approximation 
purposes. 

Our paper will be structured as follows: at Section 121 we will define our general notion 
of Young integral in the plane. Then, we will solve the wave equation at Section IHl 
Section inm is devoted to the extension of the Young integral to the light cone. We show 
how to rotate the wave equation at Section IT^ and then settle our fixed point argument 
at Section IH731 The explicit application to the fractional Brownian noise is left for Section 
13.41 Eventually, Section 0] has to be understood as a justification of the rotation trick for 
our wave equation: we explore briefiy another strategy consisting in solving the equation, 
without previous rotation of the axes. This leads to a great regularity loss of the Young 
integral, as well as some too restrictive assumptions on the driving noise. 

Along the paper we will use the notation C for any positive real constant, indepen- 
dently of its value. 

2 Two-dimensional Young integrals 

This section is devoted to a general result on Young integration in the plane, which, to 
our knowledge, cannot be found in the literature, in spite of being quite elementary: we 
consider a rectangle R = [si,S2] x [^1,^2], where Si, 52,^1,^2 are arbitrary real numbers 
such that Si < S2 and ti < t2, and we show that, under some regularity assumptions on 
the functions x,y : R ^ M., the integral J Jjiy{s,t)x{ds,dt) may be defined as a Young 
integral. 

Let us be more specific now about the regularity we will impose on the functions x 
and y, and let us define the function spaces we will consider in the sequel: first of all, for 
7,7 6 (0,1), set 

C^'^ = {/GC(i?), 11/11,,^ < 00}, (4) 



where 



■"v 'V — oU-U I II I ^ 



We also suppose that the rectangle R is contained in a sufficiently large square R = 
[i?i,i?2]^, which will be fixed throughout the discussion. With these notations in mind, 
the assumptions on x and y will be the following: 

Hypothesis (H) The function x belongs to the space C"''^' and y belongs to C'^, with 
7 + p > 1 and 7 + p > 1. Moreover, there exist two positive constants K, K' such that 

\y{s,t)-y{s',t)\ <K\s-s'\'', s,s' e [si,S2], te [tuti], 

\y{s,t)-y{s,t')\<K'\t-t'f, sG [si,S2], t,t' e [h,t2], 
with a > 1 — 7 and /5 > 1 — 7. 

Let us also define the following functional spaces, in which the solutions to our equa- 
tions will live: 

Definition 2.1. For a function y satisfying conditions (H), we define the semi-norm 

\\y\\ ■= lbllp,/3+ lbl|l:Q+ lll/||2:/3, (5) 

where the last two terms in the right-hand side denote the Holder norms with respect 
to the first and second variable, respectively. Let then C^'i be the space of continuous 
functions y such that \\y\\ < +00, and observe that we will mostly consider the particular 
case W'P := C^j. 

Let us describe now the discretization procedure we will use in order to define our 
integral on i? = [si,S2] x [^1,^2]: for any rectangle Q = [s,s'] x [t,t'] and any function g 
defined on Q, the rectangular increment oi g on Q will be defined, as usual, by 

^q9 = gi-s, t) - gi-s, t') - g{s', t) + g{s\ t'). 

For all 5 > 0, we consider (n^)^ a family of partitions of the rectangle R whose meshes 
goes to zero when 5 decreases to zero. Moreover, we assume that any of the partitions 11'' 
is formed by rectangles whose sides are parallel to the plane axes. Set 11'' = ((sf, t^))jj, 
where si = s^ < s\ < • ■ ■ < s^^s = S2, ti = t^ < t{ < • ■ ■ < tj,s = t2- With these notations 
in mind, we consider the Riemann approximations 

j=0 j=0 

where we have used the notation if^^ = [sf, sf^^] x [t^,t^^;^]. 

Before stating our basic result on convergence of Riemann sums, let us give an ele- 
mentary property concerning the partitions H^: 

Lemma 2.2. Let Ri < s < t < R2 and let s < ri < ■ ■ ■ < Vk < t be a partition of (s, t). 
Then, if k > 2, there exists an integer I G {1, 2, . . . , A;} such that 

|n+i -n-i| < tK- s\, 

with the convention that vq = s and r^+i = t. 



Proof. It is an immediate consequence of Lemma 2.2 in JH]- ^ 

We are now in a position to state the main result of this section, which gives the 
convergence of the Riemann sums defined above to a hmit zr = J J^ y{s, t)x{ds, dt): 

Proposition 2.3. Recall that we have set R = [si,S2] x [ti,t2]- Then, under Hypothesis 
(H), the sequence iz^ ) converges, as 5 decreases to zero, to some limit denoted by zr. 
Furthermore, if we consider z as a function of Si,S2,ti,t2, one gets that 

y{s,t)x{ds, dt) 

R 



< C{\\y\\oo + \\ymx\UMs2 - s,r{t2 - 1,)\ (6) 

and in particular, z defines a continuous function 



Rx R — > R 

((Sl,tl),(s2,t2)) I > Zr. 



Proof . Fix 6 > and R = [si, S2] x [ti, ^2]- We will develop the proof in several steps, as 
follows. 

Step 1. We proceed, as in the proof of Proposition 2.1 in fHI, by a kind of backward 
induction on the number of points of the partition, but instead of suppressing only one 
point, we will eliminate a whole column of 11'^. Namely, owing to Lemma 12.21 we can 
choose an integer ?e{l,2,...,fc'^ — 1} such that 

(4i-4-i)<;^(s2-si). (7) 

Consider now the new partition 11 of i? defined by 

n := {{si tj), 2 = 0, 1, . . . , ? - 1, ? + 1, . . . , /c^ J = 0, 1, . . . , k'}. 

Then, if we denote by z^ the Riemann sum corresponding to the partition 11, we obtain 
that 

z^' -z^=Yy(sL^,t^^A,s x+ Yy(stt^,)Ars X- Yy(sL^,t^^Aus .s ^.,,s ,s ,x 



R = J2 yisU,f^)Ajs ^x+Y. yislf^)Ajsx- Y. yisU,f^)A^,s ^,.^ ix[t*,t^^j 

j=0 j=0 

Y{y{slt'^-y{sU.t'^)A,s^ 



k^ — l 
i=0 



In order to get some upper bounds on z^ ~ z^, let us rewrite the last term in the above 
equality as a one-dimensional Riemann sum: set 

hl{t):=y{slt)-y{sl„t), t e [t,,t2], 



Thus, with this notation, we get 

j=0 



where vr denotes the partition of [ti,t2] given hy {ti = tf^ < tf < ■ ■ ■ < t^g = t2} ■ We will 
use now the same kind of arguments as in \1^ Proposition 2.1] to get suitable bounds on 
J^ij^)- Indeed, applying Lemma EUl we can choose an integer j G {1, . . . , A;'' — 1} such 
that 

{tU~tU)<^^{t2~t,). (9) 

One is then able to construct a new partition of [^1,^2] in the following way: 

^ := {ti = 4 < ■ • • < t|_i < 4+1 < ■ ■ • < 4 = h}. 

Hence, owing to the definition of h and /, Hypothesis (H) and the bounds ((Zj) and (jHl), we 
obtain 

I J^(7r) - J\^)\ =\hl{t^) - hl{tl,)\ X \ll{tl,) - /,^(t|)| 



|Aj5 y\ X \Ajs x\ 



< 



U,|| Jlrll ~(q^-q^ Y(f^-f^ \P(J ^ J\^(f^ -f^V 



We can now proceed in a similar way to P^ Proposition 2.1]: suppressing a carefully 
chosen point of tc and reiterating the process, one obtains, thanks to Hypothesis (H), that 

Thus, going back to (jHI), we get that 

I n^ n p5| ^ ^ \\y\\p,p\\^h,i ( \p+j(+ + \p+i 

\^R - Zr- I^i\ S<^ ,j^g _ ^y_^_^ [S2 - si) [t2-ti) , 



where E^ = hf{ti) (/f(t2) — ^f(''^i))- Now, one may reiterate the process and suppress a 
suitable column of H, so that we end up with 



j=o 



fc"*-! 



r=l 

(10) 

Let us estimate the last term of the right-hand side of the above inequality: for any 
r = 1, . . . , fc'' — 1, the integer v is an element of {1, . . . , A;*^ — 1} such that (sf^+i — sf^_i) < 
^ 's2 — si). Hence, by Hypothesis (H), one easily gets that 






r=l 



Plugging this bound in (fTUjl. we obtain that 



fc*-i 



j=0 



S f-6 iX 
J' 3 + 



+ C\\y\\i.J\x\\^,^is2 - sir+"(t2 - h)\ (11) 



and we remark that the above bound (jllll is vahd for any rectangle R and any partition 
11'' satisfying the underlying conditions. 

Step 2. Define now Z^^^ ^^ := z^, where Rq = [0, s] x [0,t], for s,t G M such that Rq C R. 
We aim to show that the sequence (Zf jn)^ has a convergent subsequence. However, owing 
to Ascoli's Theorem, it is sufficient to prove that for any positive real number K, there 
exists r] > such that 

(12) 



sup \Z'p-Z'q\<K, 
\P-Q\\<v 



where P = {s,t), Q = {s',t'). 



In order to check the above condition, we consider P = (si,ti) and Q = (82,^2) 
two points satisfying si < S2 and ti < ^2, so that we can consider the rectangle R = 
[si, S2] X [^1,^2]- Let (n'')^ be a family of partitions of the rectangle [0, S2] x [0, ^2], formed 
by rectangles whose sides are parallel to the axes. In order to avoid tiresome notations, we 
denote again by {sf, t^,), i = 0,1, . . . k^, j = 0,1, . . . , k^, the elements of 11'^. Moreover, for 
the sake of simplicity we may assume that there exist two integers l<n<k^,l<r<k^ 
such that s^ = si and t^ = ti. Then, owing to (|TT|l for the rectangles R, [0, s^_]^] x [t^_,_;^, ^2] 
and [s^_,_i,S2] x [0,t^_]^], we obtain that 



\z'-z 



n-2 k"-! 



fc"-! r-2 



i=0 j=r+l i=n+l j=0 

fc* — 1 
<C{{S2- S,r+\t2 - t,Y+^ + (S2 - Si)^+"(t2 - Uf] + 



+ c{{si_,r^\t2 - i^.r^-' + {■i_,v^%h - i^.f] + 



+c{{s2- si^.r'^itur^'' + (.2 - 4+i)"+"(d)^} + 



j=r 

5: y(0,tJ)A[o,,, 



X 



L ji j + lJ 



j=r+l 
r-1 

j=0 



X 



(13) 



Let us bound now the terms in the right-hand side of ()13p: first, the sum of all the 
terms containing products of powers of {sf — sf;)'s and (t^ — t^/)'s can be bounded, up to 
constants, by 



(S2 - s^r+^h - hr+^ + (S2 - s^r+-{t2 - t.y 



(14) 



\P+7 



(S2 - Si 



17+0 



Hence, we are left with the terms in the right-hand side of ()13|) involving sums. Let us 
sketch the calculations for the first of these terms, namely 

¥-1 

j=r 

since the two remaining terms can be treated analogously. On the other hand, S^ can be 
rewritten as a Riemann sum, in the following way: 

j=r 

where l{t) = x{s2,t) — x{si,t). Then we use the suppressing point argument, as it has 
been done at Step 1, and the regularity properties of x and y. This easily yields 

l^'l < C\\yh:p\\x\\,^4s2 - s.yih - hV^" + \y{srM)\\l{t2) - l{ti)\ (15) 

< C'||?/||2:/3||2;||^,^(S2 - Si)^ {t2 - tif^^ + C||y ||oo ||2;||^,^(s2 - Si)^ {t2 - ti)"^. 

Hence, plugging (fTlj) and (fT5|) into (fT3|l . we get an upper bound of the form 

{s2-SrY + {t2-hf 



\Z'p-Z'q\<C 



with /i, ^ > 0, from which inequality (fT2|) easily follows. Thus, owing to Ascoli's Theorem, 
there exists a subsequence of {Z^)s converging uniformly to some continuous function Z on 
R. We make an abuse of notation and we denote also by {Z^)s the underlying subsequence. 

Let us go back now to the definition of z, and for R = [si,S2] x [^1,^2], set zr = 
Zp + Zq — Z(^si,t2) ~ Z(^s2,ti)- Let us show that Zp converges to zr, when 6 tends to zero. 
For this, notice that we have the following decomposition: 

us _ n^ _ ^n^ _ ^n* _ ^n* 

^R — ^Q ^P ^[0,si]x[ti,t2] ^[si,S2]x[0,h]- 

Thus, taking limit as S tends to zero in the above expression, we obtain that z^ -^ zp. 
Furthermore, by continuity of Z, we deduce that the map (P, Q) i— > zp is continuous. 

Step 3. Let us check that the limit of {z^ )s is unique. This will be proved first when 
considering a particular integrand, namely a function x defined by some rectangular 
increment of a given path y. Secondly, in the next Step 4, we will express z^ in terms of 
a Riemann sum with respect to the function x and other suitable terms, which will finally 
lead us to the uniqueness of the underlying sequence. 

We will make use of the following fact: by a slight elaboration of the calculations 
done in the preceding Step 1 and Hypothesis (L), it can be proved that, for a rectangle 
R = [si, S2] X [^1,^2], we have 



l^n -y{s^,t^)ARx\<C\\x\\^,^[\\yt,p{s2-s^y+p{t2-t^r+P 

+ \\y\\Us2 - si)^+°(t2 - hy + \\y\\i.As2 - s^y{t2 - ti)^+^) . (16) 



9 



For a given function y satisfying Hypothesis (H), define then the function x : i? — i> M as 
follows: 

X(s,t) := A[,^,^]x[ti,t]y, (s,t) G i?, (17) 

and notice that x depends on our particular choice of rectangle R. Let us first study the 
regularity properties of x'- for ^"^Y rectangle Rq contained in i?, it can be easily checked 
that A/jgX = A/JqI/. Owing to Hypothesis (H), this implies that x ^ C'^ ■ Oii the other 
hand, if s, s' G [si, S2], s < s', and t G [^1,^2], it holds that 

\x{s,t)-x{s\t)\ = \/^[s,s']^[t,,t]y\<c{s-sY{t2-uY. 

Thus, the function xi'i'^) is p— Holder continuous uniformly with respect to t and, more- 
over, one has that ||xl|i:p ^ C{t2 — ti)^. Analogously, it turns out that x(s, ■) is p— Holder 
continuous uniformly with respect to s satisfying Hxlhip < C{s2 — Si)^. Summing up, the 
function x belongs to the space T-C^'^. 

Hence, the calculations carried out in the preceding Steps 1 and 2 hold true if we 
replace y by x- ^^ particular, by (fTBj) and the above bounds for the Holder norms of x, 
we have the following estimation: 






<K{s2-s,y+^{t2-hy'-', (18) 

where we have denoted by z^^ the Riemann sum corresponding to the function x, that is 



k^-lk"-! 



j=0 j=0 



Moreover, owing to Step 2, we obtain that the corresponding sequence {Z^)s has a conver- 
gent subsequence and therefore deduce that z^j^ converges, as 6 decreases to zero, to some 
limit zr; notice that, in order to simplify notation, we do not point out the dependence 
of X in Z^ and zr. 

Let us check that the limit of {z^fj)s is unique. For this, we follow the same lines as 
in the proof of Proposition 2.1 in [^: let Z be another limit of the sequence {Z^)s and 
set zr := Zp + Zq - Zi^^uti) - ^(s2,ti) (recall that P = (si,ti) and Q = (^2,^2)). By (UHl), 
we obtain that 



\ZR - Sr\ < 2K{S2 - s^y+P{t2 - t 



\l+p 



1, 

which is indeed true for any rectangle R. Thus, for any partition Hq = {Rij = [si, Sj+i] x 
[tj, tj+i],i = 1, . . . ,k,j = 1, . . . ,1} oi R, the following relation holds true: 

k I 

\zr - 5ij| < XI E \^R^, ^ ^R^J I 

< 2/^x 5:(s,+i - ..)^+^(t,+i - t,y^^ 

i=\ j=l 
< 2K{S2 - si){t2 - h) (^sup(s,+i - s,)^+''-i) Uup(t,+i - t.y+f'-') . 

10 



Since 7 + p > 1 and 7 + p > 1, the above supreniums tend to zero as the mesh of Ho 
decreases to zero, which proves that zr = zr. Therefore, the hmit of {z^^s is unique. 

Step 4- Now we will show that, going back to the notations of Step 3, the sequence {z^ )s 
has a unique limit. Recall that 

i=0 j=0 

where x and y satisfy Hypothesis (H). Then, the key point of our strategy is to decompose 
z^ in the following straightforward way: 

j=0 i=0 

^(19) 
where the function x is defined as in (ITTjl . Now, the uniqueness of the limit of z^j^ has 
been established in the previous step. On the other hand, owing to Hypothesis (H), we 
are able to apply the convergence results in the one- dimensional setting (see [28J) in order 
to obtain that 



Ji^ Y. 2/(^^'^i)A[sf,s^+jx[ti,t2F = j^ y{u,ti)d{x{u,t2) - x{u,ti)) 



where these limits are uniquely determined as one- dimensional Young integrals. Going 
back to relation ()19p . this finishes the proof of the uniqueness of the limit for the sequence 
[z^ )s. Moreover, the following relation is fulfilled: 

y{u, v)dx{u, v) = J J {^is^MxltiMy) ^^("' ^) + / ^^-^i' ^)'^ (^(-^2, v) - x{si, v)) 

rs2 

+ / y{u, ti)d (x(m, ta) - x{u, ti)) - y(si, ti) A^jx. 

J SI 



Step 5. Eventually, owing to ()l(i|l . it is readily checked that 



R 



y{s,t)x{ds, dt) 



<C{\\yU + \\y\\)\\x\\^,^{s2-s,r'{t2-t,r 



for any rectangle R = [si,S2] x [^1,^2], which ends the proof. D 

Remark 2.4. Assume that the functions x and y satisfy the same hypothesis as in the 
previous Proposition \2. '^ and let R = [si, S2] x [^1,^2] be a rectangle. Then, as a concequence 
of Equation Hid]) , we obtain the following estimation, which will be repeatedly applied 



11 



throughout the proof of the existence and uniqueness of solutions to the wave equation 
(Theorem\3.4\): 



y{u,v)x{du, dv) 

R 



< c\\x\\^^aA \\y\\oois2 - siy{t2 - tiy 

+ \\y\\ {is2 - Sl)^+^(t2 - tl)^+^ + (S2 - Si)^+"(t2 - t,r + is2 - S^yit2 - tl)^+'^) }. 

(20) 



3 The wave equation 

Let us turn now to the equation of main interest for us, that is the following formal version 
of a perturbed wave equation 

o2 A/" o2 A/" 

^;r(^' ^) - ^^(^' ^) = ^(^("' ^))^(^' ^)' fo^ (^' ^) ^ [0' ^] >< ^' (21) 

with initial conditions given by 



dY 

Y{0,t) = ^— (0,t) =0, teR. 

OS 

Recall that we assume that the real- valued functions X, Y are defined on [0, T] x M, 
where T is a fixed positive real number, and that the coefficient a is a real- valued smooth 
function (whose exact smoothness will be specified later on). Recall also that we give a 
rigorous meaning to equation (PT|) by means of its mild formulation, as follows: we will 
say that the continuous function F is a solution to flT^ if for any {s,t) G [0,T] x R, it 
satisfies the relation 



Y{s,t)= a{Y{u,v))X{du,dv), (22) 

J JC{s,t) 

where C{s,t) denotes the open light cone with vertex (s,t) and projected to the t— axle, 
that is the triangular domain delimited by the points (s,t), (0,t -|- s) and (0,t — s) (see 
Figure^. We will also assume that the integral defining Equation (j^^ is understood in 
the Young sense given by our Proposition 12. HI 

Before going into the details of the definition of our equation, let us specify first our 
assumptions on the function a. In fact, the main property we will need on this coefficient 
can be summarised as follows: 

Hypothesis (L) o" : M ^ M is a smooth function preserving the regularity properties on 
spaces of the form Ti^'P, for p,p > 0. Moreover, it satisfies the following two conditions: 

My)\\<C\\y\\{l + \\y\\), (23) 

IWiVi) - (^{y2)\\ < K {\\yi - y2\\oo + llz/i - 1/211) 

X (l + Ill/ill + II1/2II + llz/i - I/2II + (Ibill + llz/i - I/2II)') , (24) 
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for any y, yi,y2 G TC^'^ and some positive constants K, C. 

This assumption will be made throughout the paper, and one should observe that it is 
satisfied in the following simple case: 

Lemma 3.1. Assume that a is hounded, belongs to the space C^{M.) and has bounded 
derivatives. Then, Hypothesis (L) is fulfilled. 

Proof. In this proof, we will use the same kind of arguments as in JT], and thus only the 
main ideas of our strategy will be sketched. To begin with, we show that a preserves the 
regularity on spaces of the form Ti'''^, which amounts to control all the norms used in (0) 
to define 11 af 



First of all, let us deal with the regularity of a{y) on rectangles. For any rectangle 
R = [si, S2] X [ti, ^2], the following equality is fulfilled: 

^R<y{y) = dr dTdrdra{a{r,T)), 
Jo Jo 

where 

a{r, t) = y{si, ti) + r{y{s2, ti) - y(si, ti)) + T{y{si, ta) - y(si, ti)) + rrAny. (25) 

Thus, 

^R^{y) = I dr I dT[a'{a{r,T))drdra + a"{a{r,T))dradra\ , 
Jo Jo 

and, as in [TT], one can deduce that 

\^R^{y)\ < C\\y\\{l + \\yms2 - s,Y{t2 - t.y. 



The (1 : p) and (2 : p)— Holder regularities follow from the Lipschitz property of a. Indeed 
it is straitghforward to check that ||cr(|/)||i:p < C|||/||i:p and ||cr(|/)||2:p < C||?/||2:p. Hence 
we conclude that a{y) belongs to TC'^ and, moreover, that condition (jSHI) is fulfilled. 



On the other hand, we have to check that a satisfies the local Lipschitz property ()24|1 . 
Let us sketch the calculations concerning the (1 : p)— Holder norm; the (2 : p)— Holder 
norm may be carried out using the same arguments: let s,s,t be such that (s,t), (s, t) 
belong to R. Then, following the same lines as in ^T], it can be proved that 

(riyiis,t))-a{y2is,t))-a{yi{s,t)) + a{y2{s,t))= dr dTdrdra{¥{T)), 

Jo Jo 

where &^'(t) = y''"(s) + T{y^{s) —y^'{s)) and y"^ = 2/2 + r{yi — ^2)- Expanding the right-hand 
side of the above expression and using the assumptions on a and yi,y2, one ends up with 



I|cr(yi)-cr(y2)||i:p < C\\y2\\ {\\yi - y2\\oo + \\yi -2/2||) + \\yi - y2\\ (1 + \\yi -y2\\oo) ■ (26) 

As it has been mentioned above, we have an analogous bound for the (2 : p)— Holder 
norm. 
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Eventually, in order to deal with the regularity on rectangles of <j{yi) — cr{y2), we 
notice that the following equality holds true: 

Ai?(o-(yi) - cr(y2)) = / dr dr dudrdrd^a{a''{T,u)), 

Jo Jo Jo 

with a*" = Oi + r(a2 — cti) and a, defined as in (|^ but with y replaced by yj, i = 1,2 (see 
|llj). Then, it can be shown that 

IWiyi) - (^{y2)\\p,p < c\\yi - 1/2II 

+ C{\\yi - y2\\oo + \\yi - y2\\) {\\yi\\ + \\yi - y2\\) (1 + ||yi|| + \\yi - 2/2II) • (27) 

Putting together the bounds pH|l and (jTTj) . we conclude the proof. D 

We are now ready to define rigorously our wave equation (J22|) . 



3.1 Extension of the integral to the hght cone 



A first step towards a rigorous definition of Equation (j22j) is to extend slightly our defini- 
tion of Young integral in order to cover the case of a triangular domain like C{s, t). This 
will be done by a straightforward limiting argument, as follows. 

Let 7l{s, t) be the set of families of rectangles {Rn)n of the form [s", S2] x [t", tg], such 
that ^n=iRn = C{s,t) and 



E(4-^"r(^2-^^r<+oo. (28) 

n=l 

Then, ii X,Y satisfy Hypothesis (H) from Section 1, the integral 

a{Y{s,t))X{ds,dt) 

Rn 

is well defined, for all n > 1. Moreover, by relations © and (PHj). for {Rn)n in "^(5,^), 
the series 



a(Y(s,t))X(ds,dt):=y a(Y(s,t))X(ds,dt) (29) 

cis,t) :^^J Jr„ 

is finite and the limit does not depend on the chosen element of TZ{s,t). From now on, 
the integral J JQ/g^^aiY{s,t))X{ds,dt) will be understood by means of (P^. which gives 
a reasonable definition of a Young integral on C{s, t). 

3.2 Rotation of the wave equation 

As it will be pointed out at Section ^ it is convenient to deal with the existence and 
uniqueness of solution to Equation (pTjl after a change of coordinates corresponding to a 
45° rotation, which we proceed to detail now. Set y{s,t) := Y {^,~7§) and x{s,t) : = 
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(0,f+5) 



{0,t-s) 




Figure 1: On the left, the backward hght cone with apex {s,t) is represented, while on 
the right one can see a graphical representation of a rectangular increment of the rotated 
light cone C. 



X (^§, T/f )• Then, a trivial change of variables in the integral Equation (|^ yields that 
y satisfies 



y{s,t) 



C{s,t) 



a{y{u,v))x{du, dv), 



(30) 



where now C{s,t) corresponds to the light cone with vertex (s,t) and projected to the 
line {(s, —s), s G M}, that is the triangular domain delimited by (s, t), (s, —s) and (— t, t). 
Notice that the domain of definition of our original equation was included in an arbitrary 
large rectangle R = [Ri,R2]'^, and let R be the image of R under the rotation. Then 
we assume that the norm ||a;|| of x can be controlled suitably on R, and the domain of 
definition of the rotated equation will be 



Dn 



{(s, t) eR,-s<t<-s + V2T}, 



for a given arbitrary T > 0. The new initial conditions are given by 



y{s,-s) =0—{s,-s) 



dy 

it' 



s, — s), s e 



(31) 



(32) 



Obviously, Equation (p?n|) assumes implicitely that x and y satisfy Hypothesis (H), and 
all our statements will make use of this hypothesis. Then we will show at Section 13.41 
that this assumption can be made when X is an infinite dimensional fractional Brownian 
motion. 

Remark 3.2. Suppose that the functions x,y satisfy Hypothesis (H) and that a satisfies 
Hypothesis (L). Then, owing to Proposition \2.,Si and the considerations in the preceding 
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Suhsection VJ.li the integral of y with respect to x, namely 

I(s,t) := / / a(y(u,v))x(du,dv), 

J Jc{s,t) 

is well defined. Moreover, it can be easily seen that, for any rectangle R = [si, 82] x [^1,^2]? 
we have 

ArI = / / a{y{u,v))x{du,dv); 



IR 

a graphical proof of this fact is given in Figure Taking into account relation ^, this 
implies that I G T-C^'^ . 

Example 3.3. Let {s,t) G R and consider C{s,t) the rotated light cone with vertex {s,t). 
Assume that x, y satisfy Hypothesis (H) with 7 + 7 > 1 (this condition will be implied by 
the assumptions in Theorem \3.4^ . Assume also that o satisfies Hypothesis (L). Then, the 
integral I{s,t) is constructed as in Section 3.1. Indeed, it can be easily shown that C{s,t) 
can be recovered by the family of rectangles l+l^;^Mfc, where Mk can in turn be written as 
a union o/2^^^ squares of side ^. Then, in this case we have that 

n=l k=l \ Z / 

fc=i 
and we obtain the following estimate, which will be useful in the sequel: 

_ a{y{u, v))x{du, dv) < C\\x\\,^^ (1 + |b|| (1 + \\y\\)) {t + sy+\ 

C{s,t) 

3.3 Existence and uniqueness of solution 

We are now ready to state and prove our main general result on existence and uniqueness 
of the solution for Equation (J22I), under the form (jHDI). Recall that R denotes the image 
of -R = [Ri, -R2]^ under the 45° rotation, and that Dt has been defined by (jSH)- 

Theorem 3.4. Assume that the function x belongs to C^'^^Dt), for some 7,7 G (|,1); 
and that a satisfies Hypothesis (L). Then, there exists a unique solution to Equation IIS0\} 
in TC'^''^{Dt), for all k E (1 — 7,7) «^<^ k G (1 — 7,7), denoted by {y{s,t), {s,t) G Dt}, 
satisfying the initial conditions / T^) . 

Remark 3.5. The previous theorem ensures that there exists a unique solution to Equation 
iS(]\} in the bounded domain Dt = {(s,t) E R, —s < t < —s + ^/2T}. However, we aim 
to have a solution in the whole domain, namely 

{y{s, t), {s, t) eM.'^,-s <t<-s + V2T}. 

One may construct this solution as follows: let {Rn)n be a family of squares in M^ such 
that 
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(i) Rn C Rn+i, for din, 

We denote by yn the unique solution of Equation liSO^) on {{s,t) G Rn, —s<t<—s + 
\/2T}. Then, it is readily checked that the following is fulfilled: if n> m, 

yn{s, t) = ym{s, t), for {s, t) G Rm satisfying - s <t < -s + \p2T. 

This let us define our global solution without ambiguity: y{s,t) := yn{s,t), for all (s,t) G 
Rn with —s<t<—s + \/2T. 



Proof of Theorem\3.4\ Let us fix 7, 7 > | and k E (1 — 7, 7), k G (1 — 7,7). We will make 



use of a fixed-point argument. For this, we define the map F, from T-['^''^{Dt) into itself, 
as follows: 



r(l/)(s,t):= / /_ a{y{u,v))x{du,dv). (33) 

J Jc(s,t) 

Notice that F is well defined and, thanks to Remark 13.21 and the fact that k < 7 and 
K < 7, F(y) G 7i'*''^(Dr). The strategy in order to show that F has a unique fixed point 
is quite standard: first we prove that F maps some closed ball of 1-L'^''^{Dt) into itself and 
secondly, that it is a contraction. For the latter to be fulfilled, it is sufficient to show that 

||F(yi)-F(t/2)||<i^||t/i-Z/2||, 

for all yi,y2 G H'^''^{Dt) and some positive constant K < 1. We will develop the proof in 
several steps, which may be summarised as follows: 

1. We will prove first the existence and uniqueness of a fixed point in a small part of 
the domain, namely in a band Df, for some sufficiently small f < T. 

2. In order to iterate the procedure and cover the whole domain Dt, we will consider 
a stairs domain Sf contained in Dr. Then we will proceed to prove the existence 
and uniqueness of a fixed point in one of the squares determined by two steps of 
the stairs, denoted by Q. For this, we will recover Q by a suitable finite family 
of rectangles {Rn)n and prove that there exists a ball in 7i^''^{Rn) which is left 
invariant by F, for all n. In Figure 121 the stairs domain Sf, together with the 
square Q and the sequence {Rn)n, are represented. This step will contain most of 
the technical difficulties of our proof, since Hypothesis (L) only assumes that a is 
a locally Lipschitz coefficient, which is usually considered as a too mild assumption 
in the Young integration theory. 

3. In Step 3 we proceed to show that the corresponding map defined on the invariant 
ball of 7i'^''^{Rn) is a contraction. This will be carried out again by recovering Rn 
by a convenient family of equally sized squares. 

4. Finally, we will iterate this procedure in order to get a unique fixed point of F in a 
larger stairs domain than Sf. This will let us cover the whole band D^. 
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Step 1: Fixed point in a small part of the domain: the hand D^ 

We show first that there exists t < T such that the closed ball of T-C'^''^{Dt), namely 

B^^K := {y e n^'\Dr); y{s, -s) = 0, ||y|| < K}, 

is invariant under F, for some ^ > 0. For this, recall that, under the standing assump- 
tions, 

||y|| = lblUA+ l|y||i:K+ Wvh-.k, 

and we will try to bound all the norms above separately: 

First, let us deal with the Holder norm in rectangles. If i? = [si,S2] x [ti,t2] is a 
rectangle included in Dr, then we know that 



\^R^iy)\ 



R 



a{y{u,v))x{du, dv) 



Furthermore, owing to (|T6|) and (J23|l . we have that, 
\^RT{y) - a{y{siM))^Rx\ 

< c\\x\\,^4y\\{i + \\y\\){{s2 - s,r-^it2 - tir-") (s, - s,nt, - t,r 

< C||a;||^,^|||/||(l + \\y\\)9{T){s, - s,nt, - t.f, 
where g{T) tends to zero, as r decreases to zero. Taking into account that 

\a{y{si,ti))ARx\ < \\(y\\ao\\x\\-ia{s2- siy{t2-tiy 
and K < 7, K < 7, we obtain 

||r(y)|U,.<C^(r)(l + ||yl|(l + l|y||)), 



(34) 



where we still denote by (^(r) a positive function decreasing to zero, as r tends to zero. 

Concerning the Holder norm with respect to the first variable, we observe that, for 
any s' < s and t such that (s', t), (s, t) e Dr, we have 



\T{y){s,t)-T{y){s',t)\< 



[s',s\x[-s',t\ 



(7{y{u,v))x{du, dv) 



+ 



C{s-s') 



a{y{u,v))x{du, dv) 



(35) 



On one hand, by the same calculations carried out to obtain (jH^ or, equivalently, by (PUJ) . 
one easily gets that 



ls',s]x[-s',t] 



cr{y{u, v))x{du, dv) 



<C{l + \\yUl + \\y\\))~g{r){s-sy 



(36) 
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with g converging to zero as r \ 0. On the other hand, owing to Example 13. 3| it holds 
that 

/ / a{y{u,v))x{du,dv) < C||x||^,^(l + ||i/||(l + ||i/||))(s - s')^+^ 

J Jc{s-s') 

< C(l + IblKl + \\ymV2Ty--^'^{s - sT, (37) 

and plugging (j!?Hjl and (jTTj) in (jSHI), we obtain the following estimate: 

\\r{y)h...<C{l + \\y\\{l + \\y\\)rg{T), (38) 

where we use again the same notation ^(r) for a function converging to zero as r \ 0. 
Using the same arguments as for the (1 : k)— norm, one can also get the bound 

\\Tiy)h..,<C{l + \\y\\{l + \\y\\mT), (39) 

with g satisfying again lim^\o^('i") = 0. Therefore, putting together the three bounds 
dni), dSHl) and (jnni), we end up with 

\\T{y)\\<C{l + \\yUl + \\y\\))G{T), 

where lim^\oG'(T) = 0. Then, for any large constant K, we may choose a sufficiently 
small r such that C{1 + K{1 + K))G{t) < K. Hence, we have that, for some small enough 
r < T, 

||r(y)|| < K, whenever \\y\\ < K, 



which obviously means that F maps the closed ball B^ ^ into itself. 

We show now that F satisfies a Lipschitz property on Bfj^- C 7i'^''^{Df), with some 
f < T and Lipschitz constant K < 1. Indeed, owing to the same kind of arguments as 
before, together with Hypothesis (L), it can be proved that the following estimate holds 
true: 

\\Tiy,)-Tiy2)\\<Cir,K)GiT)\\y,-y2l 

where G{t) tends to zero as r decreases to zero. Thus, there exists a sufficiently small 
f < T such that K := G{t, ^)||a;||^,;yG'(f) < 1. Moreover, we may choose f < r. 

Now, all the previous considerations allow us to conclude that F has a unique fixed 
point in H'^''^{Df). 

Step 2: The stairs domain Sr and extension to the square Q 

We begin this part of the proof by defining what we understand by the stairs domain 
Sf. First, for a given {u,v) G D^, we denote by Pr{u,v) the open subset of the line 
r := {t = —s} n R corresponding to the projection of the rotated fight cone C{u, v) on r. 

In step 1 we have proved that Equation (j^ has a unique solution in H'^''^{Df). In 
particular, one has existence and uniqueness of solution in the subdomain Sf C D^, 
described as follows: let {Pi = {si,ti),i = 1, . . . , N} be a family of points lying on the line 
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Q 



[a,b) 



R 

^1 


"u 


«„„.. 







s s 



hV2fl 



{a.b) 



Figure 2: On the left-hand side, a graphical description of the stairs domain Sr is given, 
while on the right-hand side the decomposition of the square Q in terms of the rectangles' 
sequence (-Rn)n is represented. 



{t = — s + \/2r} such that Pr{si,ti) ripr{sj,tj] 



0, if z 7^ j, and rORc U^iPr-(si, U). 
Under these conditions, it is clear that Sf := UfLiC{si,ti) forms a stairs domain (see 
Figure El). Observe that, since we are considering a finite domain R, we may choose a 
fixed finite number A^. 

We consider now a square Q determined by two consecutive cones of 5*^, say C{si,ti) 
and C(sj+i, tj+i), for some i, as it is shown in Figure El In order to simplify notations, 
we denote by (a, 6) the right-upper vertex of Q and (a, 6), (a, &) the vertices lying on 
{t = — s -(- \/2f}, that is Q = [a, a] x [b,b]. Set Li and L2 the sides determined by the 
points (a, 6) and (a, 6) and by the former and (a, 6), respectively. Notice that we already 
know the solution on Li, i = 1,2, and this solution will now play the role of the initial 
condition for the equation on the extension Q of the domain. 

We aim to extend the existence and uniqueness result to the square Q, and a first step 
in this direction, on which we will focus for the remainder of this step, is to study the 
invariance of balls in T-C^'^IQ) under F. The main idea is to decompose the square Q in 
rectangles of the form i?„ := [s„, Sn+i] x [b, b], for n > Uq, with Uq some positive integer, 
Sno = a and Sn+i — Sn = ji (see FigureEJ. Then, we will show that in each space ?i'^''^(i?„) 
there is an invariant ball for the corresponding operator F^. In the next Step 3 we will 
focus on the contraction property of the map F^ and thus deduce the existence of a unique 
fixed point. For a given n > uq, we will use the notation || ■ ||„ = || ■ ||K,K,n+ || • ||i:fi,n+|| ■ ||2:K,n 
to denote the Holder semi-norm 11 ■ 11 defined for functions on R„. 



To begin with, we focus our attention first in the small domain Rn^ G Q G Dt- In 
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this case, we are interested in showing the existence of a unique fixed point for the map 
Tno defined for regular functions y and (s, t) G Rno as: 



(rnol/)(s,t) :=0„o(s,t) + / / , a{y{u,v))x{du,dv). 

In the previous relation, the initial condition 0„g is given by 

where the first and second term on the right-hand side of the above equality correspond 
to the known solution on the segments Li and L2, respectively. The term ■?/'„(, (s„q, 6) is 
the initial condition in the corner (s„q,6) of i?„Q. Notice then that ||0noll < 00, which 
leads us to introduce the following set of functions, on which we will build our fixed point 
argument: 

n^'^iRno) = {y ■ Rno -^ K; y\L, = ^i.no^ y\L2 = V'2,no, \\y\\no < 00} • (41) 

Observe that, for sake of clarity, we have changed a little the definition of our functional 
spaces with respect to the spaces Ti'*''', by including the initial condition in the very 
definition of 7i:'^'^(i?„J. 

Let us fix now a positive number /i such that fi < k A k and /x < 7 — k. Let also d be 
a positive constant satisfying 

Uno\\no<dn^. (42) 

Our next task is now to prove that the ball Bn^ is invariant under r„g, where Bng is defined 
by: 

Bno := {y e 7i:"'"(i?„o); b - KLo < <}• 

We will thus have to study the norm ||r„(,?/ — 0no||no; and let us first consider the Holder 
norm (1 : k). Let then (s,t),(s',t) G Rn^, with s' < s. By the very definition of the 
operator F^^ and Equation (PUJ) . we have the following estimation (notice that the initial 
condition (pno cancels out in the first inequality below): 



a{y{u,v))x{du, dv) 

s',s\x[b,t\ 



I (r„ol/) (S, t) - 0no i.S,t)- {Tn^y) (s', t) + 0^0 (s', t) \ 
< 

< Cis - s'y + C\\y\\n,a + \\y\U ({s - s'r'-'^ + is- s'r) . (43) 

On the other hand, since y G Bng and we assume condition (021), "we have that 

Ibllno < IIZ/ -0no lino + Il0no lino < {d + l)n^ . (44) 

Plugging this bound into (PHj) and taking into account that s„„+i — s„,) = ^, it turns out 
that 

||r„oZ/ - 0nol|l:.,no < "^(l + (^^ + 1)<(1 + (d + l)n(;)) = 0«). (45) 



n, 
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Analogously, we obtain the following estimate for the Holder norm WTn^y — 4>no\\2:k,no'- 

llrnol/ - </>nol|2:«,no < ^(l + {d + lH{l + {d + 1)0) = o{n^) . (46) 

Eventually, let us deal with the Holder norm on rectangles. Let R = [si, S2] x [ti,t2] be a 
rectangle included in Rn^. It is readly checked that Ancpno = 0. Thus, owing to ^K^ . we 
have that 

+ C||y||„„(l + \\y\\no)Us2 - S^y+^{t2 - tif + {S2 - ^1)^(^2 " tl)^+' 



Hence, by (B^ . (jUjl and the fact that y G Bn^, we end up with 

||r„oZ/ - </'nolUAno < ^(1 + (^+ 1)<(1 + (rf+ 1)0) = oK). (47) 

Putting together (jl5|) - (J47jl . we obtain that the ball Buq is invariant under F^^ for no large 
enough. 

At this point, let us anticipate a little on the next step, and assume our contraction 
arguments have lead us to the definition of a unique solution up to the rectangle Rn such 
that i?„ n Q 7^ 0, for n > no. We will then try to use an induction argument in order 
to define an invariant ball under the map r„+i. Since s^+i — s„ = -, even if the size of 
Rn decreases, we will cover the whole square Q in a finite number of steps, thanks to the 
fact that ^n~^ is a divergent series. Observe then that, if the solution y to our equation 
has been defined up to i?„, and if L'^'^^ denotes the left vertical side of -R„+i, then the 
solution to (pl|) on Rn+i should satisfy 

Z/|Li = V'l.n+i, and y^^n+i = ip2,n+i, 

for the function ip2,n+i = l/|_R„nL"+^' ^^^ where ^i,n+i has been introduced at relation ^^. 
We will thus introduce a space H'^''^ (Rn+i) analogously to the case n = uq given at (jH]): 

H'''''{Rn+l) = [y : Rn+l -^ K; y\LinRn+i = V'l,n+1, 2/|L^+i = ^2,n+l, WvWn+l < Oo} . 

Assume now that, for any uq < m < n, the operator F^ : 'H'^''^[R^) -^ 'H'^''^{Rm), defined 
by 



{Tmy)is,t) ■= (f)mis,t)+ , a{y{u,v))x{du,dv), 

for y G 7i'^''^{Rm) and {s,t) G Rm, leave the following ball invariant: 

Bm = {ye W^\Rn^)-^ \\y - 0„||„ < n'^}, 

where 0^ is defined as in ()4U|). We aim to show that the same is true on Rn+i- This can 
be achieved using the same kind of calculations as for the case n = hq and applying the 
following result: 
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Lemma 3.6. For all n > Hq, it holds that 

ll(/>«,IL < dn^. 



Proof: Our statement holds true for n = n^, by hypothesis. Assume that we have proved, 
for some n > n^, that ||0a:|U ^ dk^ for any k < n, and let us prove the property for (pn+i- 

Firstly, it is straightforward to check that, for all k > hq, the function 0^ has null 
rectangular increments, namely ||0A:|U,K,fc = 0. Moreover, it holds that ||0„+i||i:K,n+i < 
||'0i,n+i||i:K,n+i- On the other hand, observe that the side Li does not vary when we jump 
from n to n + 1, for all n > Uq, which implies that V'l.n+i = V^i,no- Hence, the norm 
||V^i,n+i||i:K,n+i may be bounded by a constant independent of n, say Cq. 

Let us seek now some estimates for the quantity ||02,n+i||2;K,n+i- Let s G (s„+i,s„+2) 
and t,t' G {b,b). Then, by the very definition of 0„+i, we have the following estimation: 

|0n+l(s,t) -0„+i(s,t')l = \i'2,n+liSn+l,t) - ^2,^+1 (Sn+1, t') | 



< i\\yh,k,n{Sn+l - SnT + || V^2,n ||2;K,n j (t - t')", 

where Un '■= [sn, Sn+i] X [t'jt] and y denotes the unique solution to our equation in 
Ti'^''^ (Rn) ■ Therefore, making use of the induction hypothesis, we end up with 

||^2,n+l||2:K,n+l < — -— + du'' . 

Thus, we have proved that 

||0„+i||n+i < Co + — - + dn>' < din + 1)'', 

which concludes the proof of the lemma. D 

Summing up, we have defined a finite sequence of rectangles (-R„)„ such that Q C U„-R„ 
and we have proved that there exists a ball i3„ in 7i'^''^{Rn) which is invariant under r„. 
In the next Step 3 we shall address the contractivity properties of r„. Notice again that 
we have assumed in advance that the equation could be solved on any of the rectangles 
Rk for A; < n in order to define the ball Bn+i- This claim will be justified at Step 3. 

Step 3: Contraction property 

In this part of the proof, we fix i?„ any of rectangles covering the square Q and we consider 
Bn the invariant ball for r„ in 7i'*'''(i?„). We aim to prove that r„ has a unique fixed 
point in i3„. 

In order to show that r„ is a contraction, we will consider again a suitable family of 
squares {6j,j = 1,...,M) covering i?„ and having the same size. Then we will prove 
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that the operator r„, restricted to any 6'^, is a contraction, that is, r„ defined on 7i'^''^(^j) 
satisfies the following condition: 

\\Tn{yi) - r„(i/2)llM, < K\\yi - y2\\n,9„ (48) 

for all 1/1,1/2 ^ 'H'^''^{6j) and for some constant K < 1; \\ ■ ||„e. denotes the corresponding 
norm on 'H'^''^(9j). 

Recall that Rn = [s„, s„+i] x [b,b]. Let 9i be the square contained in i?„ defined by 
[sn, a] X [b, b], for some a, b (See Figure EI). In this case, we are interested in the operator 
r„0j, defined on H'^''^{9i), as follows: 



(r„,eil/)(s,t) = 0„,ei(s,t) + / / . a{y{u,v))x{du,dv), 

J J[Sn,s]x[b,t] 

(s,t) G 9i, where (l)n,ei corresponds to the initial condition. Let us prove that if 9i has a 
sufficiently small size, then r„ g^ is a contraction. 

To begin with, let us deal with the Holder norm in rectangles. Namely, let yi,y2 G 
T-C'^''^{9i) and R a rectangle contained in 6*1, so we study the following expression: 



|A/j(r„A(yi) -r„,A(y2))| 



R 



(a(?/i(u, v)) - cr{y2{u, v)))x{du, dv) 



Owing to (j^ . Hypothesis (L), the fact that ^1,2/2 G i3„ and Lemma ESI one has the 
following estimate: 

||r„,ei(2/i)-r,A(?/2)IU,Mi < \\yi-y2\\n,e,c{n,^i){d-snr''''{b-br''''- (49) 

Concerning the Holder norms (1 : k) and (2 : /?), one uses similar arguments as for 
the above norm on rectangles to end up with bounds for ||r„^e^(?/i) — Vn,e-i_{y'i)\\i--K.,di ^^^ 
||rn,ei(|/i) — r„_0^(|/2)||2:K,6'i of the Same type as (jlHI). Therefore, if the size of 6^1, say 
(5 := a — s„, is sufficiently small, then we have that condition pH|l . for j = 1, is fulfilled. 

Arguing as for the square 9i, one could see that we can progressively cover the rectangle 
Rn by a finite family of squares {9j,j = 1, . . . , M), such that each 9j has the same size 6 
and the corresponding operator Tnfi is a contraction on 7i'^''^{9j) (notice however that 6 
depends on n). Thus, this let us conclude that the map r„ : H'^''^{Rn) — > T-C'^''^{Rn) has a 
unique fixed point. 

Step 4-' Extension to the whole domain 

Putting together the considerations of Step 2 and Step 3, we have constructed now, with 
a finite number of steps, a unique solution to (jHUjl on the whole square Q. Analogously, 
we will be able to obtain the same result for the other squares determined by the stairs 
domain Sf, in such a way that we have proved the existence and uniqueness of solution to 
Equation (J30|) in an extended stairs domain S'^, for some ri > f. Eventually, we iterate 
this procedure in order to cover the whole domain Dt of definition of our Equation ()3()|1 . 
This concludes the proof of the theorem. 

D 
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3.4 Application to the fractional Brownian motion 

In this section we apply Theorem IH.4I in the particular case where the function x corre- 
sponds to the path of some random perturbation. Namely, we are interested in Gaussian 
random noises having a fractional time correlation and some spatially homogeneous one. 
Let us make this rigorous, as follows. 

Fix H G (|, 1) and consider, on a given complete probability space (r2,jF, P), a 
L^(f2)— valued centered Gaussian process {X(0),0 G P(M^)}, where ©(M^) stands for 
the space of test functions, with covariance functional given by 

E iXi(j))XN))) = ch I dudv\u — v'^^~'^ / dxdyQ(x — y)(f)(u,x)^(v,y). 

J[0,T]2 7R2 

The value of ch is H{2H — 1) and we will focus on the case in which Q is a Riesz kernel, 
that is Q{x) = tAj, for some u G (0, 1). 

This section will be decomposed as follows: first, we will define a process denoted 
by {x{s,t), {s,t) G Dt}, which will correspond to the 45° rotation of X and will be the 
driving motion of Equation (jSO)). Then, we will study the regularity properties of x, 
namely the Holder regularity in rectangles, and therefore deduce that Theorem 13.41 can 
be applied, which will give a proof of Theorem ll.il Notice that, in order to obtain the 
C^'^— regularity of x, for some 7, 7, we will apply the following extension of Kolmogorov's 
lemma (see [Zj). 

Lemma 3.7. Let z = {z{s,t), (s,t) G M^} be a two-parameter stochastic process. Assume 
that there exist p > 1 and a,h E (-, +00) such that 

\\^Rz\\Lvin)<C{s2-s^r{t^-t^)\ 

for any rectangle R = [si, S2] x [ti,t2] and some positive constant C. Then, the process z 
admits a continuous modification whose trajectories belong to the space C" ~p' "p, for all 
a' e (i,a) andb' e (^,6). 

Let us begin now with the construction of the rotation of X: we denote by TZe. the 
45° rotation on the plane, and for any function (p G V{Dt), we set 

X(^):=X(y,o7^|). 

It is straighforward to check that {X{ip),ip G V{Dt)} is a well defined L^(i7)— valued 
centered Gaussian process with covariance functional 



E [x{ipi)X{ip2)) =CH j ^ dudv 



u-v?""-' 



. x-u x+u\ y-v y+v 
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Now, we aim to define X{lfi), where R is some rectangle contained in Dt- For this, we 
consider a sequence of functions {ipn)n in V{Dt) such that ipn{s,t) —>■ lij(s,t) as n tends 
to infinity, for all (s, t) G Dt- Then, by bounded convergence, it follows that 



lim E{\X{ipn)-X{ifm)n= lim E(|X((^„ - (^„)|^) = 0, 

?i,m— >oo n,m^oo 

and we set 

X{ln) := L^n) - lim X((^0- 

It can be easily seen that this definition does not depend on the particular approximating 
sequence. Moreover, we still have the following equality, for two arbitrary rectangles 
i?i,i?2 oi Dr. 



E (X(lji )X(ljiy)) = ch I dudvlu — V 



2H-2 

UUUf |U — v\ 

[0,T]2 



xj^^dxdy\x-y\ 1^.(^7|-'^1^.(^^7|-.^- (50) 

It is then natural to define, for {s, t) G Dt, 

x{s,t) := X (l([o,s]x[o,t])nDT) > (51) 

with the convention that, if either s or t is negative, the corresponding interval will be 
[s, 0] or [t,0], respectively. Recall that x will be the driving noise of equation (jHUj) . 

Let us fix an element u of the probability space {Q, !F, P) on which the process X is 
defined. In order to simplify the notation, we shall still denote by {x{s, t), (s, t) G Dt} the 
path of the process x associated to u. Then, the following result caracterises the Holder 
regularity on rectangles of the function x, which is obviously an important step in order 
to apply Theorem 13.41 

Lemma 3.8. Almost surely on fi, the function x : Dt — > M defined by A51]} belongs to the 
space C"^'^ , with rj G (0, 7) and fj G (0, 7), for any 7, 7 such that 7 + 7 = i7 + ^y^. 

Proof We will apply Lemma IT7I For this, we fix a rectangle R = [si,S2] x [ti,t2] in 
Dt and we compute the square moment of A/jx. Indeed, by the very definition of x, we 
obtain that 

E{\Anx\') = E{\Xiln)\') = E{\X{ln o R.)]')- 

Hence, owing to (jKUjl . we end up with 



£■ MAijxpj =ch / dudv\u — V 



2H-2 

auavya — v\ 
[0,T]2 



y. I «.|. - yrU [^. -^) 1. (^, ^j . (52) 

and notice that, for instance, we have 1/? (~^5 ^^) = '^ni^^^)^ where R := TZ-^{R). 
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In the sequel, we will make use of the following notation: 

and set also As := S2 — Si and At := ^2 ~ ^i- With these notations in hand, we will try 
to get a bound of the type 

E (lA^xp) < C{^sY\^tf\ (53) 



where 7, 7 satisfy 7 + 7 = // + ^y^, for a general rectangle i?. However, throughout the 
proof, we will assume that As < At: indeed the case At < As may be treated analogously, 
and when i? is a square, relation (j3^ follows easily. 

In order to get good bounds of the right-hand side of (jK^ . we will decompose the 
indicator function 1^, as follows (see Figure Ej): 

where 

• Ri is the triangle determined by 

{(5*11, Til), (5'i2, T12), (5*12, 5*12 - y2si)}, 

• i?2 is the parallelogram determined by 

{(5'i2, T12), (S'12, 5*12 - V2si), (5*21, T21), (5*21, 5*21 - V2S2)}, 

• -R3 is the triangle determined by 

{(5*21, T21), (521, 5'21 - y2S2), (5*22, T22)}. 

Owing to this decomposition, the integral in the right-hand side of ()52|) shall be split in 
six terms. 

The proof can be divided in two different situations: 
First case: At < 3As. Under this assumption, it is straightforward to check that, in the 
terms of the integral ()52|1 involving the indicator function of the parallelogram R2 -which 
are the most winding-, there will always be intersection between the domains of variation 
of X and y. Moreover, the fact that As < At < 3As will let us obtain the appropiate 
bounds for ()52j) without too much effort. The details for this case are omitted here. 
Second case: At > 3As. Here, fixed u and v in their respective domains of variation, 
there does not need to be intersection between the domains of variation of x and y in the 
terms involving the parallelogram i?2- As it has been mentioned, it turns out that this 
situation leads to much more complicated calculations in comparison with the above one. 
Therefore, we will just deal with the proof in this second case. 

As mentioned before, our partition of R into three pieces leads to the computation of 
six different terms. However, by symmetry, it is sufficient to consider the following three 
situations: 
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(Sn-^n) 



i^u'^n-^^ ) 



(^21. ^n) 




{5,,,5,,-V^.J 



(5,,.^,,) 



(Si.'^iJ 



Figure 3: Decomposition of the domain R = TZ-il{R). 



Term corresponding to lj:^^{u,x)lj^^{v,y) 



We have to get bounds, for instance, of the integral 

dy{u — V 



du dv dx -"-'- -^2^-2' 

•512 •! Sl2 JU—y/2s2 JV—y/2s-i 



x-y\ 



(54) 



the case when v G (u, 5*21) can be treated in the same way. The main idea is to decompose 
the domain of integration above so as to get rid of the absolute value there, and therefore 
be able to compute, or eventually bound, the remaining integrals. Notice that, for the 
sequel of this proof, we have decided not to write the multiplicative constants which 
appear in our estimations, and which can vary from line to line. 

First, it holds that (jH^ equals to 

/•5i2+\/2As fu fu—\pls\ pv—\/2s\ 

Ai + A2:= du dv dx dy{u - vf^''^\x - y\''' 

J S\2 '^ S12 Ju—-\/2s2 Jv—\f2s2 

pS2l pu pu~y2s\ pv—y2s\ 

+ du dv dx dy{u - vf^~'^\x - y\~'' . 

JSi2+V2As JSi2 Ju-V2s2 Jv-V2s2 

Then, on one hand we have that 



Ai= du dv dx dy{u-vY^'^{x-yy 

J S12 J S12 J u~\/2s2 Jv—\f2s2 

(■Si2+\/2As pu pv—\/2si pv—\/2si 

+ du dv dx dy{u-vf^-'^{y-xy 

J S12 J S12 Ju—\/2s2 J X 

/•5i2+\/2As fu /■u—\/2si /•v—\/2si 

+ du dv dx dy{u — vY^~'^{x — y)~^ , 

J S12 J S12 Jv—\/2s\ Jv—\f2s2 
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and it is readily checked, from this decomposition, that 
On the other hand, one has that A2 = Bi + B2, with 

rS2i /•u—\/2As /•u—\/2si /•v—\/2si 



El := du I dv dx dy{u — v) {x — y) " ^ 

B2 := / du I dv I dx I dy{u — v) ^ |x — y]^''- 

Notice that, in the domain of integration of Bi, the following estimation holds true: 

(x - yy <{u-v- V2As)-''. 

Hence, 

Bi < (As)' / du / dv{u - ^;)2^-^ 

"'Sl2+V^As JSl2 

and computing the above integral, one ends up with the bound 
Eventually, it turns out that we can decompose B2 as follows: 

rS2i ru /•v—\/2si rx 

B2 = du dv dx dy{u — vY^~~'^{x — y)~^ 

^Si2+v^As Ju-V^As Ju-\/2s2 Jv-\/2s2 

/•S'21 ru rv-\/2si rv-\/2si nu n 

+ du dv dx dy{u - vf^'^{y - x)'" 

JSi2+V^As Ju-\/2As Jm-V2s2 Jx 

/■S'21 /■" /■«— v2si pv—\/2s-i 

+ du dv dx dy{u — vY^~'^{x — y)~'^ , 

JSi2+V^As Ju-V^As Jv-V^si Jv-^/2s2 

and from this expression, it is not difficult to prove that 

B2< iAsf^-''+\At + As). 

Putting together all the bounds that we have obtained on Bi, B2 A2 and Ai, we conclude 
that the term (jH^ can be bounded by 

^/^sYH-u+2 ^ (As)2(As + At)''"-'' + (As)2^-'^+^(Ai + As). (55) 

Let 7, 7 belong to (|, 1) and such that 7 + 7 = if + ^y^. Then, owing to (jKKjl . one can 
easily see that (jH^ may be bounded, up to constants, by (As)''^(At)''^. 

Term corresponding to lpf_^{u,x)lff^{v,y) 
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Let us treat now the term corresponding to lj^_^{u,x)lj^,^{v,y), namely 

du dv dx dy{v — u) ~ \x — y\~'^. (56) 

Su JSi2 J-u+V2ti Jv-V2s2 

This term equals to Di + D2, where 

rSi2 fSi2+V2As fu-V2si fv~V2si 

dy{v — u) \x — y\ '^, 

S12 rS2i /■u—\/2si /•V—V2S1 



D2 = du dv dx dy{v — u) {y — x) " ^ 

JSu JSi2+\/2As J~u+\/2ti Jv-V2s2 

and using similar arguments as for the term Ai, one can easily show that the following 
estimate holds true (recall that we omitt to write the multiplicative constants): 

D2 < {As)\At)^^-\ 

For the term Di, we can show, along the same lines as for the term A2 above, that 

Di < (As)2^-'^+2_ 

Thus, dSni) niay be estimated by 



27/'Af^27 



(AsYHAt) 
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with 7 + 7 = H + ^. 



Term corresponding to Ijj (n, a;)lj^ {v,y) 



Our aim now is to bound the integral 

'•S12 rS\2 j-u-V2si j-v-V2si 

du I dv 1 dx 1 dy\u — v\ ^ \x — y['^ . 

5ii JSii J-u+V2ti J-v+V2ti 

It turns out that this term involves easier calculations in comparison with the two other 
situations, and that it can be bounded by 

{AsY^-'^^ < {AsY^AtY\ 

The details of computations for this term are left to the reader. 

Putting together the estimations obtained in the study of the three terms above, we 
get that 

E{\Anx\')<CiAsY^iAtY\ 

for a positive constant C. Taking into account that we are in a Gaussian context and 
applying Lemma ITTj we get the desired regularity for the process x. D 

Remark 3.9. As mentioned before, Theorem M.lM s now an easy consequence of the previ- 
ous lemma. Indeed, it allows us to apply Theorem \3.4\ to the fractional Brownian motion x 
defined by / I5i)) . which in turn implies Theorem M . IV by the simple rotation trick of Section 
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4 Sharpness of our method 

Let us come back for a moment to Section \'A.2\ and to the definition ()H()|1 we gave for the 
wave equation: we performed a rotation TZz. on our initial equation in order to get a new 

domain of integration C{s, t) whose sides are parallel to the axes. Then, as pointed out in 
Remark 13 .21 if x, y : Dt —>■ ffi are two functions satisfying Hypothesis (H) and a : M ^ M 
fulfills Hypothesis (L), the integral 



I(s,t) = / / a(y(u,v))x(du,dv) 

J JC(s.t) 



IC(s,t) 

is well defined in the Young sense. Moreover, it has the same Holder regularity on rect- 
angles as the function x. 

However, as it will be made clearer later on, one could define directly the integral 

I(s,t) =11 a{Y{u,v))X{du,dv) = T f Gs-uit,v)a(Y(u,v))X(du,dv), 

J Jc{s,t) Jo JR 

where the fundamental solution G of the wave equation is given by 

Gs{t,v) = -l{it-v\<s}, s>0,t,v eR, 

by means of Young approximations with respect to a good partition of the domain, such 
as a dyadic one. In this section, we will try to show that this direct strategy does not 
behave as well as the one we proposed at Sectional Indeed, we will show that, even in 
the linear case (i.e. a = 1), we are quite far away from obtaining for X the same Holder 
regularity as the control X. Let us make this rigorous, as follows. 

Let us first consider the linear case, that is we want to define the Young integral 



Iis,t)= G,^uit,v)X{du,dv). (57) 

Jo Jm. 

We have the following result on existence and regularity of the integral X. 

Proposition 4.1. Recall that C^'^ has been defined by relation ^, and assume that X 
belongs to the space C''^ , with 7 + 7 > 1. Then, the expression \51\) is well defined in 
the Young sense. Moreover, the function {X(s,t), (s,t) G [0, T] x M} belongs to C''', with 
77 = P7 + 7 - 1 andfj = {1- p)^, for all p G (^, l) . 



Remark 4.2. According to the preceding result, we deduce that 77 + ^7 = 7 + 7 — 1. Thus, 
the integral! has substantially lost regularity with respect to the control X . Furthermore, 
here and in the remainder of the section, the Young integrals we consider are not really 
defined as such, since they are based on a particular partition, suitable for computations. 
We have chosen that solution for sake of simplicity, but we believe our calculations could 
be carried out for a general family of partitions too. 

Let us sketch the proof of the above proposition. 
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Proof of Proposition \4 1\ We use a Young approximation to define the integral in the 
right-hand side of (J57jl . Fix (s,t) G [0, T] x M and consider the following dyadic grid on 
the rectangle in M^ determined by {(0, t — s), (s, t — s),{0,t + s),{s,t + s)}: 



u7 = s^. i = OA 2", 



/, J-, . . . , ^ , 



^; = ^-^ + si^> J = 0,1,...,T 



Notice that the above rectangle contains the support of the function {s,y) t-^ Gs^u(t,v) 
and each of the rectangles of the grid is of size ^ x ^. We set also 



and denote by Gs,t{u,v) the function Gs-u(t,v)l[Qs-^{u). 

The natural Riemann type sum approximating the integral in (jFfj) would be 

Equivalently, we study the convergence of the series 

oo 

J2iJn+l{s,t)-J„{s,t)), (5^ 

n=l 

which we will decompose in the following way: first, notice that 

A/" X = Arn + l X -\- Arri + l X -\- Arll + l X + AfJl + l X. 

«J 2i,2j ^2i+l,2j ^2i.2j + l -'2i+l,2j + l 

On the other hand, the term J„+i(s, t) can be written as 

2"-12"+l-l r 

Thus 

Jn+i(s, t) - J„(s, t) = Ai{s, t, n) + ^2(5, t, n) + ^3(5, t, n), 

where 

2"-l 2"+l-l 

A(.,t,n) = 5: 5: (a,K+\<^)-G,,K+\,t;27'))A,,^.^X, 

2"-l 2"+i-l 



A,{s, t,n) = Y: E (Gs.K+S <^) - a,tK^\ <+\)) ^q^.^^x. 

i=0 j=0 '' -^ 



i=0 j=0 

Next, using the fact that Gg^t is the indicator function of the light cone, a possible strategy 
for the estimation of those terms is the following: 
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1. Identify the non-trivial contributions in the sums Ak{s,t,n), k = 1,2,3, and then 
bound them in a suitable way. 

2. In order to get the desired bounds for the non-trivial contributions, use the fact that 

- -^ \ m(7+7) 



IArXI < C , 

for a rectangle R of order m in the grid defined above. 
Using these ideas, one can obtain the following estimation: 

C 



|J„+i(s,t) - J„(s,t)| < 



2"(7+7-l) ' 



with some positive constant C independent of n. Since by hypothesis 7 + 7 > 1, we 
deduce that the series (J58|l is convergent, and its sum defines properly the integral in the 
right-hand side of (jHTjl . 

Let us outline the main computations leading to the regularity result for the function 

{X(s,t), (s,t) e [0,r] X M}: we fix < Si < S2 < T and -00 < ti < ^2 < +00, and we 
will study the quantity |ArX|, where R = [si,S2] x [ti,t2]- Furthermore, we assume that 
ti + S2 < ^2 — S2, which implies that the light cones with vertices (s2, ti) and (s2, ^2) do not 
intersect (the other case can be treated analogously). With these assumptions in mind, 
we have 

A/jJ = J(S2, ^2) - ^{S2, h) - I{si,t2) + J(si, ti) 

or, equivalently, 

/ / {Gs2,t2{u,v) - Gs2,tAu,v) - Gs^^t2{u,v) + Gsi^ti{u,v)) X{ds,dy). (59) 

Jo JR 

We set 

D{u,v) =Gs2,t2{u,v) -Gs2,ti{u,v) -Gsi^t2{u,v) + Gs^^tAu,v), 

and we proceed to the analysis of ()59j) along the same lines as in the first part of the 
proof. Namely, we will use Young type approximations of the following form: fix Uq G N, 
and let io G Z_ be such that 



C:=^i + (^2-.i)^<0 and <Vi:=^i + (^2-.i)^ 



We consider the dyadic grid on [0,T] x R given by 11" = 11" l+l 112, where 11" is the grid 
on the rectangle in M^ determined by {{u'l°,ti - S2 + u^°), (s2, k- S2 + u'^°), {u^°,ti + S2 - 
Ui°), {s2,ti + S2 - u^°)}, I = I, 2, in the following way: 



< = si + (s2-si)-, 2 = 2"-"%,..., 2", 



^f = ti + {S2 - si)^, J = -2"(1 - 2-"«zo), . . . , 0, . . . , 2"(1 - 2-"o^o) 
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for all n > riQ. Observe that each square of the grid 11" is of size ^^;r^ x ^^;r^, / = 1, 2. 
Set now 



jl,n 



Ln Ln 



in ra 1 ., r 'i" 'i" 1 



and notice that, once again, we will try to get some information about the convergence 
of the series 



7 , {Jn+l{Sl, S2', ti, t2) — Jn{Sl, 82] ti, ^2)); 



(60) 



n=no 



where 



Jn{Sl, S2', tl, t 



Y: D{u^,v^pA^,„X. 



Now, the desired convergence can be obtained in the same spirit as for the first part of 
the proof, but with slightly more effort. For sake of conciseness, the details are left to 
the reader, but let us briefly justify the need of introducing the parameter p. Namely, it 
turns out that one of the terms that we have to deal with for the convergence of (jHUIl is 
given by 

2 —In — I / 

Ain) = y (A.2,n+1 X-A.l,n+1 X 

,_. V 2i„ + l,2fe -'2i„ + l,2fe 

where in is some integer such that is contained in (■u"^,M2j^^i]. In order to get bounds 
on A{n) of the form (s2 — Si)^{t2 — ti)'^, p, ft E (0, 1), we notice that 



A ,2,n+l X — A J 



2i„+l,2fc 



where 



Ri 



l,n + l 
2i„+l,2fc 



„ n+1 n+1 

"2i„+l5 "2i„+2 



X — /\f/nX — Ar" X 



-^fc+l/2 ' 



X 



l,n+l 2, n+1 



Then, owing to the hypothesis on X and the definition of (m", f /"), for any p G (0, 1), we 



obtain 



2"-i„-l 

\Mn)\< E 



Ar2,n+1 X — Arl,n + 1 X 
2in+l,2k 2i,i+l,2fc 



/\fjnX — Apn X 

"-k "■k + l/2 



1-p 



<C2 
<C 



S2 - Si \^('>'+^) /'S2-Si\ (^-P'^'y 



2n+l 
]_\ "(P7+7-1) 

2 



(t2 - tl) 



{1-P)7 



{S2 - SiY''+^t2 - h)'''-''^\ 

This bound is then used to obtain the desired regularity in our statement. 



n 



Let us jump to the non-linear case, which amounts to give a rigorous meaning and 
obtain regularity properties of the integral 



Xz{s,t) 



Gs-u(t, v)Z{u, v)X{du, dv), 



for a given process Z. More specifically, we consider the following hypothesis for the 
function Z: 
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(H^) The function Z : [0, T] x M ^ M is bounded, Z(0,t) = 0, for all t G M, and there 
exist a G (0, 1) and a positive constant C such that 

|Z(si,t)-Z(s2,t)|<C|Si-S2r, 

for any si,S2 G [0,T], t G M. 

Then we get the following existence and regularity result: 

Proposition 4.3. Assume that X belongs to the space C^'^ , with 7 + 7 > |. We suppose 
also that the function Z belongs to Cp^i+^i-W^-pYi ^ mth p e (^ 3-27-7 ^ 27+^7-2 ^^ ^^^ ^^^^ 
hypothesis (Ha) is satisfied for alia G (2(2 — 7 — 7), !)■ Then, the integral I zis,t) is well 
defined in the Young sense and it defines a function belonging to C^^'^"'^^'^^^''^^ . 

Remark 4.4. Here again, for sake of conciseness, we will omit the proof of this propo- 
sition, which is much more involved than the one of Proposition \4 1[ but may be carried 
out using analogous arguments. 

Remark 4.5. The motivation to construct the pathwise integral Tz is to solve, given some 
control X , the integral equation 



Y{s,t)= / / Gs^uit,v)cr{Y{u,v))X{du,dv). 

Jo JR 

To prove existence and uniqueness of solution to the above equation, one usually applies 
a fixed point argument or, equivalently, a Picard iteration scheme. For this reason, in the 
preceding Proposition \4.3\ we have assumed that the function Z has the same regularity as 
the integral in the linear case (see Proposition ^. 1\} . 

Remark 4.6. According to the statement of Proposition |7^ we deduce that not only the 
integral Xz looses again regularity with respect to the control X, but also the regularity of 
the latter must be strengthened, since we assume now 7 + 7 > 5/3 instead 0/7 + 7 > 1. 
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